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Q . Abstract By studying various rational integrable maps on C with p invari- 

qq \ ants, we show that periodic points form an invariant variety of dimension > p 

for each period, in contrast to the case of nonintegrable maps in which they are 
isolated. We prove the theorem: 'If there is an invariant variety of periodic 
points of some period, there is no set of isolated periodic points of other period 
^ . in the map. ' 

O 

2 : 1 Introduction 

When a map is given there is no immediate way to foresee the fate of its iterations or 
whether it is integrable or nonintegrable. If the map is integrable the effect of a small 
difference of initial values remains small, while it behaves chaotically for some initial values 
if the map is nonintegrable. In spite of this large difference one cannot easily distinguish 
the two cases. 

It would be desirable to know a way to distinguish integrable maps from nonintegrable 
ones just by investigation of the first few steps. This article is a contribution to finding 
^ | such a criterion in the case of some higher dimensional maps. 

The notion of integrability of a map has been studied in various contexts, such as 
the inverse method, the Baklund transformation method, the Darboux transformation 
method, the direct method, the geometrical method, etc. They judge the integrability of 
a map by testing if the map belongs to the classes to which the methods are applicable. 
Since the procedure is rather difficult, they will not decide if the information is limited 
only to the first few steps of the map. 

In the case of continuous time Hamilton systems the Liouville theorem guarantees 
integrability of the systems if there are sufficient number of involutive invariants. There 
is no such criterion if the time evolution is discrete. In fact, when d — 1 invariants are 
given arbitrary, we can always derive a map which reduces, for example, to the logistic 
map after the elimination of d — 1 variables by using the invariants. 

In the theory of ordinary 2nd order differential equations, nonexistence of movable 
branch point singularities, dependent on initial values, indicates integrability of the equa- 
tion (the Painleve test). The method of singularity confinement was proposed [TJ [2] to 
replace this test and apply to discrete maps, but it does not always work [3]. As far as 
one dimensional maps are concerned there have been intensive studies mainly from the 
mathematical point of view [U [5j [6j [7] . It is known, for example, that the appearance of a 
Julia set, the closure of the set of repulsive periodic points, characterizes nonintegrability 
of the map. On the other hand very little is known about higher dimensional maps. In 



practice the best we can do at present is to study statistical information of the maps [3l [8] . 

We consider a rational map on C d , where C = {C, oo}, 

x = (x 1 ,x 2 ,... > x d ) -> X = (X 1 ,X 2 ,...,X d )=:X^. (1) 

Throughout this paper we count the degree of freedom d only which are coupled and are 
dependent each other. We are interested in the behaviour of the sequence: x — ► X^ 1 ) — > 
X( 2 ) — > • • ■ . In particular we pay attention, in this paper, to the behaviour of periodic 
points of rational maps. If the map is nonintegrable we may find a set of isolated points 
with fractal structure as a higher dimensional counterpart of the Julia set [9]. We would 
like to know what object appears when the map is integrable. 

We study in §2 the general feature of the periodicity conditions = x of a rational 
map with p (> 0) invariants. After the elimination of p variables there are d — p such 
conditions for each period. If they are independent, we obtain a set of isolated periodic 
points in general. We call the periodicity conditions of this generic type 'uncorrelated'. 
It may happen that the d — p conditions are correlated each other. In this case we find 
varieties of periodic points instead of isolated points. If all of the d — p conditions are 
correlated, i.e., when the conditions are 'fully correlated', the varieties will be determined 
only by the invariants. We call such a variety 'an invariant variety of periodic points'. 
Their precise meaning will be explained in §2. We prove there the following theorem: 

Theorem 

// there is an invariant variety of periodic points of some period, there is no set of 
isolated periodic points of other period in the map. 

Although the theorem doesn't tell us directly whether the map is integrable or nonin- 
tegrable, it is remarkable that it discriminates clearly a particular type of maps with an 
invariant variety of periodic points from the rest. In generic case with no invariant variety 
the periodicity conditions are uncorrelated and the periodic points form, for each period, 
a set of isolated points. In fact, if a rational map is given arbitrary and we calculate the 
periodicity conditions of any period, we find them uncorrelated almost always. 

Now we recall that almost all nonlinear maps which are chosen arbitrary are noninte- 
grable in general. This means that the chance for a nonintegrable map not to have a set 
of uncorrelated periodicity conditions of any period must be practically zero. A typical 
higher dimensional nonintegrable map has a higher dimensional analog of the Julia set. In 
this case the periodic points form a closure of fractal sets of isolated points which become 
the source of chaotic orbits. 

The precise notion of nonintegrability of a map has not been known. This makes diffi- 
cult to discriminate integrable maps from nonintegrable ones unambiguously and prevents 
us to proceed discussion further. In order to overcome this problem we are going to adopt 
the facts we discussed above as our working hypothesis: 

// a map is nonintegrable, there is a set of uncorrelated periodicity conditions of some 
period. 



Note that we do not require that the periodicity conditions of all periods of a non- 
integrable map are uncorrelated. We also emphasize that the existence of uncorrected 
periodicity conditions of some period does not mean nonintegrability. Once we adopt the 
hypothesis, the following statement will be justified directly from our theorem. 

Conjecture: 

If there is an invariant variety of periodic points of some period, the map is integrable. 

Note that our conjecture does not exclude the possibility that some integrable maps do 
not have an invariant variety of periodic points. For example an integrable map with no 
invariant does not have an invariant variety. On the other hand there are d dimensional 
maps which have d — 1 invariants but reduce to some nonintegrable maps after elimina- 
tion of d — 1 variables. Therefore the situation is quite different from the continuous time 
Hamiltonian flow, whose integrability is guaranteed by the Liouville theorem if there are 
sufficient number of invariants. 

We are not going to prove above hypothesis in this paper. But it is not difficult to 
convince ourselves of the fact as we calculate the periodicity conditions of any period of a 
rational map which is chosen arbitrary. In order to support our conjecture we would like 
to discuss in §3 and §4 if a set of periodic points of any period forms an invariant variety 
when the map is known integrable and there are sufficient number of invariants. As far 
as we have examined, we have found that all periodicity conditions possess this property 
if the map is integrable and periodic points exist, while no such property has been found 
otherwise. We show in §3 some examples of integrable maps, such as the Lotka-Volterra 
maps and the Toda maps. We also study the g-Painleve maps which are integrable but 
not volume preserving in general. We find invariant varieties only when the parameters 
are restricted to have invariants. 

In §4 we discuss one dimensional maps to which some of higher dimensional integrable 
maps, such as the symmetric QRT map, the 3d Lotka-Volterra map, the g-Painleve IV 
map with certain coefficients, can be reduced by using invariants. Analysis of these maps 
enables us to derive iteratively infinite series of invariant varieties of all periods. 

The method developed in §4 enables us to study properties of many higher dimensional 
maps all together simply by studying a single one dimensional map. We apply this 
method in §5 to investigate the integrable-nonintegrable transition of higher dimensional 
maps. In particular we consider a simple one dimensional map to which many higher 
dimensional maps can be reduced and interpolates between the Mobius map and the 
logistic map. By studying the behaviour of periodic points in detail we clarify how the 
Julia set collapses and the invariant varieties of periodic points are created when the 
parameter of interpolation approaches to the critical value of the transition. This example 
shows explicitly the phenomenon which takes place at the border between integrable and 
nonintegrable regimes of higher dimensional maps. 



2 Nature of the Periodicity Conditions 



We study in this section the nature of the periodicity conditions of a rational map on C d 
which has p (> 0) invariants. The periodic points of period n will be found by solving 

xf = x 3 , j = l,2,...,d. (2) 

If ifi(x), i?2(x), • • • , H p (x) are the invariants, the solutions of ([2]) are constrained on an 
algebraic variety of dimension d — p specified by the set of equations 

Hi(x) = hi, z = l,2,...,p. (3) 

Here hi, hi, ■ h p are the values of the invariants determined by the initial point of the 
map. Let us denote this variety by V{h), i.e., 

V(h) = {x| # 4 (x) =ht,i = 1, 2, ...,p}. (4) 

The problem of finding periodic points is equivalent to finding an ideal generated 
by the set of d + p functions {X- n ^(x) — Xj, i/j(x) — hi}. Since the existence of the 
invariants enables us to eliminate p components of x from (12"]) , the ideal reduces to the 
pth elimination ideal generated by certain functions r^ n ) satisfying 

^{h, h 2 , h p , 6,6, -, 6- P ) = 0, a = 1,2,. ..,d-p, n>2. (5) 

Here by 6,6, •••,6- P we denote the variables which parameterize the variety V(h) after 
the elimination of the p components of x. Note that the fixed point conditions (n = 1) 
are excluded in (jSJ) since they are nothing to do with the invariants. 

For an arbitrary set of values of hi, hi, h p , the functions T^'(h,^) define an affine 
variety, which we denote by V^((r)), i.e., 

K (B) «r» = {e|rW(/»,0 = 0, a = l,2,...,d-p}, n>2. 

In general this variety consists of a finite number of isolated points on V{h), hence zero 
dimension, corresponding to the solutions to the d — p algebraic equations (jSJ) for the d—p 
variables 6,6, ■■■,£,d- P - Once the values of 6,6, •••,6-p are decided by solving (jSJ), the 
location of a periodic point on C d will be determined from the information of the values 
of hi, hi, hp. In this case we obtain a number of isolated periodic points of period n, 
and say that the periodicity conditions <^ are 'uncorrected'. Needless to say this case 
includes a map with no invariants. 



There are possibilities that the equations (jSJ) impose relations on hi, hi, h p instead 
of fixing all 6' s - Let I be the number of such equations. We write them as 

^ (h, hi,..., h p ) = 0, a = 1,2,..., I, (6) 

instead of to emphasize independence from 6's. If m is the number of the rest of 
the equations 

r^(/ il ,/ i2 ,...,/ ip ,6,6,-,6- P ) = 0, « = l,2,...,m (7) 

d — p — in variables are not determined from the periodicity conditions. This means that 
y( n )((r)) forms a subvariety of dimension d — p — m of VVh). We say that the periodicity 



conditions are 'correlated' in this case. Since the conditions leave p — I invariants free, 
there remain totally d — m — I variables undetermined by the periodicity conditions. In 
other words the points of period n satisfying the conditions (jHJ) and ([IJ) form a subvariety 
of dimension d — m — I in C d . 

When m = the periodicity conditions determine none of the variables ^1,^2, 
but impose I relations among the invariants. Then the affine variety V^((r)) coincides 
with V(h). In other words every point on V(h) is a periodic point of period n, while V(h) 
itself is constrained by the relations among the invariants. We say that the periodicity 
conditions are 'fully correlated' in this particular case. If we replace hi by -ffj(x) in 7^(/i) 
the periodicity conditions 

7i n) (ffi(x), # 2 (x), flp(x)) = 0, a = 1, 2, I. (8) 

enable us to consider the constraints on the invariants as constraints on the variables x. We 
denote by v^ n '((j)) the affine variety generated by the functions 7( ^(ifi(x), H 2 (x), if p (x)), 
and distinguish it from V^((T)). Namely we define 

7 i n) (#i(x), ff 2 (x), ff„(x)) = 0, a = 1, 2, /}. 



The significance of defining v ^((7)) lies in the fact that if a point belongs to v^ n \(j}) 
the point is a periodic point of period n. This is true only in the case m — 0. If m > 0, 
the constraints on the invariants are not sufficient to determine the periodic points. We 
also notice that the variety v ^ n \(l)) is determined by the invariants hi, h 2 , h p alone as 
is clear from the construction. We have called ^^((7)) 'an invariant variety of periodic 
points' in §1. Let us summarize the properties of f^ n '((7)) as follows. 

• The dimension of i^ n '(( 7 )) is d — I (> p). 

• Every point on ^^((7)) can be an initial point of the periodic map of period n. 

• All images of the periodic map starting from a point of ^^((7)) remain on it. 

• i>^(( 7 )) is determined by the invariants of the map alone. 

Having introduced some notions about the periodicity conditions we are now going to 
establish our theorem presented in §1. For this purpose let us suppose that the periodicity 
conditions of period n are given by ([H]) and (j7|). They are fully correlated if m — 0, and 
correlated if d — p > m > 0. Note that, since the I relations among the invariants are 
assumed independent, the conditions © are consistent only if 

p> I, d — p > I + m (9) 

are satisfied. We also emphasize that the conditions (EJ) include neither fixed points nor 
points whose periods are the divisors of n. 

We further assume that the periodicity conditions of period n) are uncorrelated, 
i.e., 

ri k \h l ,h 2 ,...,h p ,£ 1 ,&,...,Z d - P )=0, a = 1,2,..., d-p, k^n (10) 



«W«7» = {x 



and ask if there exist £1, £2, £,d- P which satisfy ([6]), ((Tj) and ({TO]) simultaneously for some 
values of the invariants h\, h 2 , h p . It is necessary, but not sufficeint, for this to happen 
that the total number d of the variables £1, £2, and the invariants h\, ...,h p exceeds 

the total number d — p + / + m of the equations (jSJ), (13) and (TTOl) . i.e., p > I + m must 
hold. Combining this together with (Q we obtain 

min{p, d — p} > I + m. (11) 

Now we suppose that (TTTT) is satisfied and ffTOl and (j7]) are solved to determine all 
£'s and m invariants, say hi, fo 2 > •••) ^m- Since the rest of the invariants, h m+ i, h p , are 
completely free it is aways possible to fix them such that (jSJ) is satisfied. Therefore it 
is inevitable, as far as I 7^ and the periodicity conditions of period k is uncorrelated, 
there exist some points which have two different periods n and k simultaneously. This is 
certainly a contradiction. Thus we are lead to the following conclusion: 

Lemma 

A set of correlated periodicity conditions satisfying ffl\) and a set of uncorrelated pe- 
riodicity conditions of a different period do not exist in one map simultaneously. 

In particular when the periodicity conditions of period n are fully correlated, i.e., 
when m = in (JTj), the condition (iTTj) is always satisfied as long as ([6]) has solutions. Our 
theorem in §1 follows from this fact immediately. 

Some remarks are in order: 

1. The existence of an invariant variety of periodic points is a sufficient, but not nec- 
essary, condition for the integrablity. 

2. The existence of a set of isolated periodic points is a necessary, but not sufficient, 
condition for the nonintegrability. 

3. When p < d/2, a set of correlated periodicity conditions satisfying d — p > I + m > 
p > I > and a set of uncorrelated conditions of other period can exist together in 
one map. 

4. When p > d/2, uncorrelated periodicity conditions of two different periods can 
exist together, although the total number of the variables exceeds the number of 
the conditions. This is because the argument, which lead us to the Lemma, does not 
apply when I = 0. Therefore there are possibilities that all periodicity conditions are 
uncorrelated even though the map has invariants more than d/2. A simple example 
will be discussed in §4. 

3 Invariant Varieties of Periodic Points 

In order to see the correspondence between the existence of an invariant variety of periodic 
points and integrability of the map, we are going to study in this section various integrable 
maps and see if the periodicity conditions (j2J) exhibit the invariant varieties of periodic 
points. 



Since we need assistance of computer, the procedure of finding the invariant varieties 
of many variables and/or of higher degree becomes harder as the dimension of the map 
and/or the degree of period increases. Nevertheless, as far as we are able to calculate, 
we are convinced that all perodic points, if they exist, satisfy fully correlated periodicity 
conditions and form an invariant variety for each period, as we present below. 



3.1 Lotka-Volterra maps 

Although there have been intensive studies of discrete integrable systems pllj HU [121 E3 
[H], the precise notion of integrability of a map has not been known. Among others the 
bilinear method provides a powerful tool to find integrable maps, in the sense that the 
solutions are given explicitly in terms of the r function of the KP hierarchy [TOl [T5l [16]. 
We study, in this section, the Lotka-Volterra maps, the Toda map, and the g-Painleve 
maps, all derived by this method. 

The Lotka-Volterra map x — > X of dimension d can be found [15j [16] by solving 

X 3 {1 - X^ x ) = Xj (l - x j+1 ), j = 1, 2, d (12) 

for X = (Xi,X 2 , ...,Xd) under the conditions Xj + d = Xj (j = 1,2, ...,d). We show in the 
Appendix that the invariants of this map are given by 

= ^j-L,j2,...,j k x h x h ' ' ' x ik(X ~ •fy'i— ~ x h~i) ■ ' ' (1 — ^ifc-i) 

(A; = 1,2,..., [d/2]) (13) 

r = x x x 2 ■ ■ ■ x d 

Here the prime in the summation J2' °f (Tl~3|) means that the summation must be taken over 
all possible combinations ji,j2, ■■■■,3k but excluding direct neighbours. The total number 
of the invariants is p = [d/2] + 1, where [d/2] = d/2 if d is even and [d/2] = [d— l)/2 if d 
is odd. 



3.1.1 3d and 4d Lotka-Volterra maps 

When d = 3 the explicit form of the map is obtained as 

l-X 2 +X 2 X 3 I-X3 + X3X1 l-x 1 +x 1 x 2 , v 

Xi = xi- ■ , X 2 = x 2 - ■ , X 3 = x 3 - ■ . (14) 

1 — X3 + X 3 Xi 1 — Xi + X\X 2 1 — x 2 + X 2 X 3 

by solving (TT2"]) for (X 1 , X 2 , X 3 ). There are two invariants 

r = X!X 2 x 3 , s = (1 - xi)(l - x 2 )(l -x 3 ), (15) 

where we used s := H\fr instead of H\ for convenience. Since d — p = 1, we expect one 
invariant variety of periodic points for each period. After some manipulation we find the 
following results for the variety 

^)(( 7 )) = {x| S + l = 0} (16) 

for the period 2 case, 

v {3 \(>y}) = {x| r 2 + s 2 - rs + r + s + 1 = 0} (17) 



for the period 3 case, 

t/ 4 )(( 7 )) = { x | r 3 s + s 3 - 3rs 2 + 6r 2 s + 3rs - r 3 + s = 0} (18) 

for the period 4 case and 

t> (5) (( 7 )) = { x | rV - rV - 6rV + 10r 3 s 6 + 3s 5 r + s 6 + s 5 + 3rV - 3rV 
-6rV - r 6 s 3 + 3rV + s 4 + 21sV + 6s 4 r + rV + s 7 + 27s 5 r 2 
-3s 6 r - rV + 21r 2 s 6 - 10r 3 s 3 - 6rs 7 + s 8 = 0}, (19) 

for the period 5 case, etc.. 

The 4 dimensional Lotka-Volterra map is defined by 

1 - X 2 - X 3 + X 2 X 3 + X3X4 1 - 2 3 - £4 + X3X4 + X4X1 

At — Xi- ■ ■ , X 2 - x 2 - ■ ■ , 

1 — £3 — X4 + X3X4 + X4X1 1 — X4 — Xi + X4X1 + X1X2 

1 — X4 — Xi + X4X1 + XiX 2 1 — Xi — X 2 + XiX 2 + X 2 X 3 

X% — X 3 , X 4 — X4 . 

1 — Xi — X 2 + XiX 2 + X 2 X 3 1 - X 2 - X 3 + X 2 X 3 + X 3 Xi 

There are three invariants of the map, given by 

r = XiX2X 3 X4, Hi = Xi + X2 + X 3 + X4 — X1X4 — X2X1 — X 3 X2 — X4X3, 

H 2 = XiX 3 + X2X4 — X 2 X 3 X4 — X3X4X1 — X4X1X2 — XiX 2 X 3 + 2XiX 2 X 3 X4- 

We find an invariant variety of periodic points 

t/ 2 )(( 7 » = {x| ^-2 = 0} 

for the period 2 case, and 

w (3) (( 7 )) = {x| Hi - H\H 2 + 2HiH 2 - 2rH 2 - 6H 2 + Hf - 2H\ 
+2rHi + 2Hi + r 2 + 2r + 5 = 0} 

for the period 3 case, etc.. 

3.1.2 5d Lotka-Volterra map 

When d = 5 there are three invariants, thus we have d — p = 2. From our argument of 
§2 the periodicity conditions ([2]) will be satisfied by imposing two independent conditions 
7 j n ) = 0, 7 2™^ = 0, simultaneously, for each period. The situation changes from the 
previous cases d — 3,4 where d — p is 1. 



Xi = xi 

X 2 = x 2 

X 3 = x 3 

X4 = X4 

X 5 = x 5 



g 




in the 


case d = 5 we obtain 






1 


- x 2 


-x 3 - 


X4 + X 2 X 3 + X2X4 + X 3 X4 + X4X5 


- X2X3X4 


- X2X4X5 + X2X3X4X5 


1 


- x 3 


— X4 — 


X 5 + X3X4 + X3X5 + X4X5 + X5X1 


- X3X4X5 


- X3X5X1 + X3X4X5X1 


1 


- x 3 


— X4 — 


X5 + X 3 X4 + X3X5 + X4X5 + X 5 Xl 


- X 3 X4X 5 


- X3X5X1 + X3X4X5X1 


1 


— X4 


-x 5 - 


Xl + X4X5 + X4X1 + X 5 Xl + XlX 2 


— X4X5X1 


— X4X1X2 + X4X5X1X2 


1 


— X4 


-x 5 - 


Xl + X4X5 + X4X1 + X5X1 + X1X2 


— X4X5X1 


— X4X1X2 + X4X5X1X2 


1 


- x$ 


-xi- 


X 2 + X5X1 + X 5 X 2 + X1X2 + X2X3 


- X5X1X2 


- X5X2X3 + X5X1X2X3 


1 


- x 5 


-xi- 


X 2 + X5X1 + X 5 X 2 + XlX 2 + x 2 x 3 


- X 5 X ± X 2 


- X5X2X3 + X5X1X2X3 


1 


- Xl 


- x 2 - 


X 3 + XlX 2 + X1X3 + X 2 X 3 + X3X4 


- XlX 2 X 3 


— X1X3X4 + X1X2X3X4 


1 


- Xl 


- x 2 - 


X 3 + XlX 2 + X1X3 + X 2 X 3 + X3X4 


~ XlX 2 X 3 


— X1X3X4 + X1X2X3X4 



T, J- Tr*Tn 4- ToTj -1- ToTj -I- T^Trr — Tn'TiTA — MiTc -I- ToT^TiTt 



From our general formula the three invariants are 



Hi = Xix 2 + x 2 x 3 + £3X4 + X4X5 + x 5 Xi — Xi — x 2 — x 3 — £4 — x 5 , 

H 2 = X1X3 + X2X4 + X 3 X 5 + £4X1 + X 5 X 2 - X1X2X3 - X 2 X 3 X 4: - X 3 X 4 X 5 

—£4X5X1 — X5X1X2 — X1X2X4 — X1X3X4 — X1X3X5 — X2X3X5 — X2X4X5 

+X 2 X 3 X4X 5 + X3X4X5X1 + X4X5X1X2 + X 5 XiX 2 X 3 + X1X2X3X4, (20) 

r = 21X2X3X4X5. 

If we form two particular combinations H2 + 3Hi + 5 and Hi + r + 2 from these invariants, 
we see that the Grobner basis of this ideal generates the 3rd elimination ideal of the 
functions {Xj — Xj}. Therefore the invariant variety of periodic points of period 2 is 

t/ 2) (( 7 )) = {x| H 2 + 3Hi + 5 = 0, Hi + r + 2 = 0}. 

This is an algebraic variety of dimension 3. 

3.2 3 point Toda map 

The 3 point Toda map [151 HBj 

(ii,i 2 ,i 3 ,Vi, v 2 , v 3 ) 

is defined by 

T . hvi + i 3 h + v 3 vi 

«2^3 + «2«3 + V 2 V 3 
. %iV 2 + %l%2 + ViV 2 

h = ^3 , 

I3V1 + t 3 ii + v 3 Vi 

T . i2V 3 + % 2 % 3 + V 2 V 3 

h — ^ l 5 

llV 2 + %l%2 + ViV 2 

This map has four invariants, 

h = i\ + i 2 + i 3 + vi + v 2 + v 3, 

t2 = i\h + «2«3 + hk + viv 2 + v 2 v 3 + v 3 vi + i\v 2 + i 2 v 3 + z 3 t>i, 
t 3 = hhh, 

t' 3 = fi^3- 

From these data we find the Grobner bases generating the 4th elimination ideal of 
{xj 2 -* — Xj} and {X- 3 -* — Xj}. The periodicity conditions for the period 2 case are not fully 
correlated but only correlated. Therefore we have no invariant variety of periodic points. 
In the case of period 3 the conditions are fully correlated. We obtain an invariant variety 
of periodic points of dimension 4 as follows, 

</ 3 )« 7 }) = {x| ti = 0, t 2 = 0}. 



{h,I 2 ,I 3 ,V u V 2 ,V 3 ) 



t r l2V 3 + i 2 % 3 + v 2 v 3 
V\ = vi . 

l 3 Vi + % 3 %i + V 3 Vi ■ 

t t- hvi + i 3 ii + v 3 vi 
V 2 = v 2 - — . 

l\V 2 + %i% 2 + flf2 ' 

t t kv 2 + %\i 2 + viv 2 
V 3 = v 3 . 

«2^3 + + V 2 V 3 



3.3 g-Painleve maps 



We examine in this subsection discrete analogues of the Painleve equations which do not 
belong to the Lotka-Volterra series. Since they do not have invariants in generic cases 
their integrability is not obvious at all. As we choose the parameters in the map properly 
there appear invariants and our theorem guarantees integrability of the map. Therefore 
they provide examples of the maps in the border to which our theorem can apply to 
discriminate integrable maps from nonintegrable ones. 

There have been proposed various types of the discrete Painleve equations which 
preserve integrability [TTl [IB"] . Among others we study here the symmetric versions. By 
choosing the dependent variables properly the g-Painleve IV map can be written in the 
symmetric form [18] 

1 - a 3 x 3 + a^xsXi 

1 — OL\X\ + 0>\Ot2X\X2 

1 - aixi + aia 2 xix 2 , 01X 

A 2 = a 2 a 3 x 3 , (21) 

1 - a 2 x 2 + a 2 a 3 x 2 x 3 

1 - a 2 x 2 + a 2 a 3 x 2 x 3 

X 3 = a 3 a t xi- ■ . 

1 — a 3 x 3 + a 3 a\X 3 X\ 

This map is integrable in the sense that the map admits Backhand transformations. It 
is also supported by the existence of some type of explicit solutions. Since the Jacobian 
of this map equals ai«2«3, the map does not preserve the volume in general. Under these 
circumstances the periodicity conditions are uncorrelated for all periods. Hence this map 
is an example of integrable maps which do not have an invariant variety of periodic points. 

As we impose the condition oio 2 a3 = 1, the map has one invariant 

r = (22) 

hence p = 1. Because p < d/2 in this case, correlated and uncorrelated periodicity 
conditions of different periods can exist in one map, as we noticed in Remark (3) of §2. 
If we solve the periodicity condisions of period 2, we find an invariant variety 

i;( a )(( 7 » = {x|r+l = 0} (23) 



and also a curve in C 3 

{x |(1 - 03^3)^2 - a 2 l {^ ~ x l) x 2 - (x 3 - a 3 )x 3 = 

n (0:10:2 - x 3 )(xix 2 - x 3 ) - oia;i(l - x\) = }. (24) 

The first solution v ^ of (|23|) corresponds to the fully correlated periodicity conditions 
of the case I — l,m — 0. According to our conjecture the existence of guarantees 
integrability of the map. The second solution, the curve of (1241 . fixes isolated points on 
V(h), corresponding to the uncorrelated conditions of the case I = 0, m = 2, as we fix the 
value of r. 

The result of period 2 provides us useful information to explore the implication of our 
theorem. The theorem tells us that the existence of excludes uncorrelated periodicity 
conditions of other period in this map. Thus (124")) is the only possible set of periodic 
points which can generate isolated periodic points for each value of r. On the other hand 
the existence of the solution (1241 excludes invariant varieties of other periods. Therefore 



the periodicity conditions of all other periods must produce surfaces in C 3 corresponding 
to I = 0, m = 1. In fact we find the points of period 3 only on the surface defined by 

o 2 (l + o 2 + a^a\){r 2 + 1) + a±(l + a\ + 3o20 2 )(x 2 X3 + rxi) 
+o 2 o 2 (o 2 + a\a\ + 3)(x 3 Xi + rx 2 ) + Oi0 2 (l + a\a\ + 3o 2 )(xiX 2 + rx 3 ) 
— (1 + o 2 + a^i) fai (0:2^1X2 + x 2 x 3 + Oi0 2 x 3 Xi) + (1 + o 2 + o 2 o 2 ))r 
— «i(l + a-L + a 2 Q; iX a 'i + tti«2^2 + ^2^3) 

1 2/2, 222, 2 2\, 2/222, 22, 2 2 2 2\ 

+a 1 {x 1 + a l a 2 x 2 + o 2 x 3 ) + a 1 (a 2 x 1 X2 + x 2 x 3 + a 1 a 2 x 1 x 3 ) 

— 2o 2 o 2 (x 2 x 2 + a.\x\x-i + oio 2 xix 2 + oix 2 X3 + oio 2 xix 2 + x 2 x 2 ) = 0. 

For each value of the invariant r this surface generates a curve on V(h). Similarly we 
have found a surface formed by points of period 4, but more complicated. 



If we further fix the parameters to 01 = o 2 = 03 = 1, the set of equations (I2T]) 
is obtained from the 3d Lotka-Volterra map (fUj) simply by the shift (Xi,X2,Xs) — > 
(X 3 ,Xi,X 2 ). This map has two invariants, one is r of (|22|) and the other one is s = 
(1 — Xi)(l — x 2 )(l — x 3 ), the same as (fT5|) . In this particular case we have found a series 
of the invariant varieties of periodic points 

w (2) (( 7 )) = {x|r + l = 0} 

w (3) (( 7 )) = {x| r 2 + s 2 -rs + r + s + l = 0} (25) 
w (4) (( 7 )) = { x | s 3 r + r 3 + 6rs 2 + 3rs - s 3 + r - 3r 2 s = 0}, 

and so on. 

The symmetric version of the g-Painleve V map has been given explicitly [19] by 

1 — 03X3 + 0:30:4X3X4 — 030:40:1X3X4X1 
; > 





= Oi0 2 X 2 


x 2 


= o 2 o 3 x 3 


x 3 


= 030:4X4 


x. 


= 0401X1 



X LXlO^l \JL\*~X 2 Uj\JU2 l-JL\K-X2*- JL 3 Uj l ,Aj Z' jU d 

1 — O4X4 + O4O1X4X1 — 040i0 2 X4XiX 2 

1 — 2 X 2 + 2 3 X 2 X 3 — 2 030 4 X 2 X3X 4 ' 
1 — Of 1X1 + Oi0 2 XiX 2 — Oi0 2 03XiX 2 X3 

1 — 03X3 + O3O4X3X4 — O3O4O1X3X4X1 ' 
1 — 2 X 2 + 2 3 X 2 X 3 — 2 3 04X 2 X 3 X4 



(26) 



r = XiX 2 X3X4, 

and an invariant variety of periodic points 

v ( 2 )(( 7 )) = {x| r-1 = } (27) 

in the period 2 case. If we further fix the parameters as 01 = o 2 = 03 = 04 = 1 there are 
two more invariants 

s = (1 - xi)(l - x 2 )(l - x 3 )(l - x 4 ), t = (1 - x 2 x 4 )(l - X1X3). (28) 

in addition to r, and (|27|) is the invariant variety of periodic points of period 2 again. 

As the period increases the computer manipulation becomes more difficult to derive 
the invariant varieties explicitly. For some cases we discussed in this section, the method 
we develop in §4 will enable us to find iteratively the invariant varieties of all periods if 
we know one of them. 



4 Series of Invariant Varieties of Periodic Points 



We study in this section a method of generating an infinite series of invariant varieties of 
periodic points. Although application of this method is limited to some particular type of 
maps at present the existence of such method seems to suggest that the invariant varieties 
of different periods are correlated each other. 



4.1 Generation of higher dimensional maps 

Let us suppose that the p invariants i?i(x), H 2 (x), H p (x) are given and we want to 
derive a d dimensional map which reduces to the following d — p dimensional one 

Xi -> Xi = fi(x!,X2, ...,x d -. p ), i = l,2,..., d-p, (29) 

after the elimination of the variables Xd-p+i, ■ ■■,x d . This will be done easily. If 

Xi(x) = gi(x), i = d-p+l,...,d 

are the solutions of 



Hi(fi,f 2 , ...,fd- p ,X d -p+i, ...,X d ) = Hi(x), i = 1,2, ...,p, 

then 

x — >■ X = (/i, f d _ p , g d - P +i, g d ) 

is the d dimensional map which satisfy our requirements. In this way we can derive many 
higher dimensional maps, irrespective whether they are integrable or not, which have a 
number of invariants. 



As an illustration, let 

1 + bx 

H(x,y) = y , b, c: constants (30) 

1 + cx 

be the invariant and try to find a two dimensional map which reduces to 

x ^ X = h^^-(l + cx). (31) 
1 + bx y ' K ' 

where h is the value of H{x,y). Following to our prescription we immediately find the 
map 

( \ (vv\ (f > \ (l + bx)(l + cy(x + a))\ 
(x,y) (X,Y) = [(x + a)y, yf— — . . 32 

\ (1 + cx){l + by{x + a)) J 

The map ( 1311) is integrable in the limit c — > 0, since it becomes the Mobius map. When 
c is finite the map ( 13~TT) is nonintegrable since it is equivalent to the logistic map. Therefore 
the map (|32l) provides an example of nonintegrable two dimensional map which has one 
invariant. In general the existence of a large number of invariants does not guarantee 
integrability. In fact we could derive a map which reduces to the logistic map, if we 
choose f\(x) = hx(l — x) in (1291) when there are given d — 1 invariants. 



4.2 Iteration of the Mobius map 



In the previous subsection we learned that there are many higher dimensional integrable 
maps which reduce to a common integrable lower dimensional map. We study in this 
subsection the periodicity conditions for the higher dimensional maps which reduce to 
the Mobius map. 

Now suppose that a d dimensional map of x is reduced to the Mobius map 

x -> X = h — , (33) 

1 + bx 



after the elimination of d — 1 variables. If the map is the one of (1321) with c = 0, the 
parameter h is the invariant y(l + bx), whereas a, b are some constants. Generally a, b 
and h are functions of the invariants, but not dependent on x otherwise. 

The iteration of this map does not change the form of the map but only changes the 
parameters a, b, h. If we write 

X in) = tfn) X + » W (34) 

1 + bWx 

after n steps, the (n + l)th parameters are related to the nth ones by 

a(n+1) = a + a^h^ b(n+1) = fr (w) + bh^ h(n+1) = h feW + a6W 

h(n) + a0 (n) ' I + ft/j,(n) a (n) ' 1 + lfa{n) a {n) ' 

from which we can determine all parameters iteratively as functions of the initial values 
(a, b, h). 

The periodicity conditions of period n for the map ( 1331) are now satisfied if the param- 
eters (a, b, h) satisfy 

(a (n+1) , 6 (n+1) , /i (n+1) ) = (a, b, h). (35) 

From our argument it is clear that the periodicity conditions do not fix the values of the 
variable x but impose some constranits on the parameters, hence on the invariants. 
Solving (|35|) iteratively we find the invariant varieties of periodic points as follows 



(2) 


= {x 


(3) 


-{■ 


(4) 


= {x 


(5) 


-{■ 


(6) 


= {x 



l + /i = 0} 

1-h 3 



1 - h 

, U2 



abh = 



h 5 

+ abh(3 + (A + ab)h + 3h 2 ) 



h 



(36) 



4.3 Biquadratic maps 

We have seen that there exist infinitely many d dimensional integrable maps which reduce 
to the Mobius map (|33l) . We notice that the Mobius map (1331) can be written in a bilinear 
form bXx + X — hx — ah = 0. If we study some of well known integrable maps with 



d — 1 invariants, however, many of them do not reduce to the Mobius map but to the one 
specified by the equation 

S(X,x;q) = (37) 

where 

S(X, x; q) := aX 2 x 2 + b(X + x)Xx + c(X - xf + dXx + e(X + x) + /, (38) 

and 

q = (a, b, c, d, e, /) G C 6 (39) 

are parameters dependent of the invariants. We shall call a map x — > X of this form a 
'biquadratic map' in the following. Some examples which reduce to the biquadratic map 
( 157)) are the symmetric version of the QRT map [HI H2] , the 3d Lotka- Volterra map of 
( TT4"j) and the g-Painleve IV map ( 12T1) with «i = «2 = «3 = 1, as we discuss later. 

One might wonder that the function (138)) does not determine the image of the map 
uniquely. Moreover the number of images increases rapidly as we repeat the map. This, 
however, does not happen because of the symmetry of the map under the exchange of x 
and X. In fact we see that one of the solutions of S(X, x; q) = is X = X^ and another 
is X = X*- -1 -* corresponding to the forward and the backward map of x. As we repeat 
the map we obtain the biquadratic map every time, but with different parameters, as we 
explain now. 

The problem of finding the image Q(x) of the second iteration of the map ()37j) is 
equivalent to finding an elimination ideal generated by S(Q,X; q) and S(X,x; q). After 
the elimination of X the function, say W^ 2 \Q, x), which generates the 1st elimination 
ideal, becomes quartic in both Q and x. Thus we obtain four solutions of Q for given 
x. Those we expect are the solutions of S(Q, X^ ±lS) ; q) = 0, corresponding to the paths 
x -> X (±1) X (±2) and x X (±1) — > x, since the X, which was eliminated, could be 
either X« or X^ 1 ). 

In fact, after some manipulation, we obtain 

W^(Q,x) = (Q-x) 2 S(Q,x;q 2 ), (40) 

with the new parameters q*- 2 -* = (a^ 2 \ b^ 2 \ c^ 2 \ d^, e^ 2 \ f^) being given by 

a (2) : = (ae - cbf - {ad - 2ac - b 2 ) (be - cd + 2c 2 ) , 

b {2) := (ae- cb)(2af - be + cd - 4c 2 ) - (ad - 2ac - b 2 )(bf - ce), 

c ( 2 ) : = ( a f -c 2 ) 2 -(ae-bc)(bf - ce), (41) 

d^ ■= A(af -c 2 ) 2 - 2(ae - be) (bf - ce) - (be - cd + 2c 2 ) 2 

-(ad - 2ac - b 2 )(df - 2c/ - e 2 ), 

e ( 2 ) : = (fb-ce)(2af -be + cd-4c 2 )-(fd-2fc-e 2 )(ea-cb), 

/( 2 ) : = (fb-ce) 2 -(fd-2fc-e 2 )(be-cd + 2c 2 ). 

Now we suppose that the kth images X^ ±k ^ of the map ( 1371) are found up to k = n by 
solving S(Q, x; q^ k ') = 0, k = 2, 3, n for Q. The (n + l)th image X^ n+1 ^ of x must be 
obtained by solving the equation W^ n+1 \Q, x) = which is derived from S(Q, X; q(™)) = 
and S*(X, x;q) = after the elimination of X. Now by solving S(Q, X^ 1 ); q( n )) = for 
Q we obtain X^ n±1 \x) and X^~ n±1 )(x). We know already X( ±(n ~ 1 ))(a;) as the solutions of 



S(Q, x; q*™ -1 )) = 0. Therefore W^ n+1 \Q, x) must have S(Q, x; q^ n_1 ^) as a factor. Since 
W^ n+1 '(Q,x) is a biquartic function and S(Q,x; q( n_1 )) is a quadratic function of Q and 
x, both are symmetric under the exchange of the variables, another factor of W^ n+l \Q, x) 
must be a symmetric and biquadratic function of Q and x. Such a function can be written 
uniquely in the form of fl38l) itself. Therefore we have 

W {n+1 \Q,x) = SiQ^-q^SiQ^;^^). (42) 

The new coefficients q( n+1 ) will be read off by comparing both sides of fj42l as functions 
of q and q^ n \ Using the notation (g A g') n = 99 — g'g^ n \ the results are as follows: 



= -^((a A c)l-(a A b) n (b A c) n ), 



b {n+1) = (w _ 1) I (w _ 1) ((a A b) n (b A c) n - (a A c) 2 ) + (a A c) n ((Q A e) n + 2(fr A c) w ) 



- ]^({ a /\ h )n{b A e) n - (a A &) n (c A d) n + (a A cf) n (6 A c), 
c ("+i) = _L_ ^( ce (n) _ 6/(n))(ae(™) _ & c (")) + ( c ftW _ ea (n) )(/6 (n) - ec (n) ) 

+ (a/ (n) - cc (n) ) 2 + (/a (n) - cc (n) ) 2 ^) , 

d (n+l) = ___ ( _ y(n-l) a (n+l) _ a (n-l) /( n+l) _ 46 («-l) e (n+l) _ 4e (n-l) 6 (n+l) + 

d^ n 1 > \ 

+ (c A d) 2 n - (a A b) n (e A f) n - (6 A c) n (c A e) n + (a A d) n (d A f) n + 2(b A e) ri {a A f) n 
3(c A e) n - (b A f) n - (d A e)nj(3(b A c) n - (a A e) n - (b A d) n ^j 

+2((a A d) n - (a A c) n )((c A /) n - (d A /) n ) 
+2((6 A c)„ + (a A e) n ) ((6 A /) n + (c A e) r 

i /p(n-l) 

= ^(j^((/Ae) n (e A c^ 

~ ^((/Ae) n (e A 6)n - U\ e )n(c A d) n + (/ A d) n (e A c) n )j, 
/ (n+1) = ^Ty((/Ac) 2 -(/ A e) n (e A c)„). (43) 



e 



Our problem of studying the behaviour of the map (|37|) has been converted to studying 
the behaviour of the parameters , & f^ n '). These parameters determine a 

trajectory of iteration of the map (l3"Tj) : 

. . . <_ X (- 2 ) <— X^ 1 ) <— x — X« — X^ 

Namely we can consider (|43p as a map q( n ) — > q( n+1 ) in C 6 . 

The periodicity conditions of period n for the map x — > X defined by (|3"T|) are satisfied 
if the parameters satisfy 

q (n+1) (q)=q. (44) 



Despite the complicated expression of the relation (1431) . we observe a special dependence 
on the nth parameters q^ n \ Besides c^ n+1 \ the dependence of the (n + l)th parameters on 
the nth ones is always in the form (gAg')n = gg'^ ' — 9 ' 9^ '• They all vanish simultaneously 
when the periodicity conditions (144p are 'fully correlated'. In other words if there exists 
a function 7^ n+1 )(q) such that 

q(»)(q) = q + 7 (n+1) (q)q (n) (q) (45) 

so that (jUJ) is satisfied by a single condition 7 n (q) = 0. Here we use the same notation 
which we used already as a function of h in higher dimension. Although they arise 
in different contexts they turn out to be the same object as we will see later. 
When (T4"5l) holds, the equation S(Q,x; q( n+1 )) = can be written as 

c ("+l) (Q _ %) 2 + 7 (n+l) ( q )^(«+l) (Q, x ) = . (46) 

Here 

K (n+1) (Q,x) = a {n+l) Q 2 x 2 + 6 (n+1) (Q + x)Qx + d {n+l) Qx + e (n+1) (Q + x) + f {n+l \ 

and cS n+l \ for instance, is obtained from a^ n+1 ^ simply replacing (g^g') n by (#A#')n- If <5 
is a point of period n + 1, the first term of (|46j) vanishes. Hence the periodicity condi- 
tion requires for the second term to vanish. This is certainly satisfied for arbitrary x if 
7 < - n+1 )(q) = 0, namely when the periodicity conditions for the parameters q*™- 1 are fully 
correlated. The other possible solutions obtained by solving K^ n+1 '(x,x) = will not 
correspond to the points of period n + 1, but represent the fixed points or the points of 
periods which divide n + 1. 

Let us present the functions explicitly in the cases of small number of n. We 

can show 

(a A b)2 = (af — eb — 3c 2 + cd) {2a? 'e — abd + b 3 ) 
(a A c) 2 = (af — eb — 3c 2 + cd) (a 2 / + ac 2 — acd + b 2 c) 
(b A c)2 = (af — eb — 3c 2 + cd){2ace — abf — be 2 ) 

(e A /) 2 = (af-eb-3c 2 + cd)(edf-e 3 -2bf 2 ), 
from which we find 7 < - 3 ' ) (q). We can derive other results iteratively as follows: 
7 ( 3 )(q) = a f - be - 3c 2 + cd, 

7 {4) (q) = 2acf - adf + b 2 f + ae 2 - 2c 3 + c 2 d - 2bce, 

7 (5) (q) = a 3 f+(-cf 2 d + 2cfe 2 + fde 2 -3ebf 2 -e 4 -c 2 f 2 )a 2 

+ ( - 13c 4 / + 18c 3 fd + de 3 b + 2cf 2 b 2 + 7dc 2 e 2 - ce 2 d 2 - 2ce 3 b 
+ 2c 2 feb - 7fd 2 c 2 - Uc 3 e 2 + cd 3 f + fb 2 e 2 + f 2 db 2 - ebd 2 f)a 

- cd 2 b 2 f - 6 3 e 3 - Ac 3 deb + cdb 2 e 2 + 13ec% - f 2 b 4 + 7fb 2 c 2 d 

+ c 4 d 2 - 5c 5 d + 5c 6 - 2fb 3 ec - e 2 c 2 b 2 + eb 3 df - Ufb 2 c 3 , (47) 

and so on. The formula (|43l) enables us to continue finding a series of 7*- n ^(q) systemat- 
ically. From the way of this construction it is apparent that all 7's are functions of the 
invariants alone if one of them is so, a fact being consistent with our theorem. 



4.4 Maps reduced to the biquadratic one 



Having established the scheme of generating the invariant varieties of periodic points of 
the map (13T1) , we are going to present some examples of higher dimensional maps which 
reduce to the biquadratic one by using invariants. If {Xj = Fj(x), j = 1,2, ...,d} is the 
map and H^ipc), i — 1,2, d—1 are the invariants, we calculate the 2(d— l)th elimination 
ideal of the set of the functions 

{Xj - i^x), iZi(x) - h h j = 1, 2, d, i = l, 2, d - 1} 

by eliminating the 2(d — 1) variables (x 2 , x 3 , x^, X 2 , X 3 , X d ). 

4.4.1 The symmetric QRT map 

Consider the two dimensional map 



v'(y)p"(y) 


- p'(y)v"(y) 


-x(p'(y)<j>"(y) 


- 0W(v)) 




- <P'(y)p"(y) 


- x(cf)'(y)ri"(y) 


- vW(y)) 



(x,y) -> (X,Y) = | [ — ) ' , — J . ( is 

Here 

4>'(x) := a'x 2 + b'x + c', <t>"{x) '■= a " x2 + b"x + c" , 

7]'{x) := b'x 2 + (d' - 2c')x + e', T]"{x) := b"x 2 + (d" - 2c")x + e", 

p'{x) := c'x 2 + e'x + f, p"(x) := c"x 2 + e"x + f", 

and q' = (a', b', c', d! , e', /') and q" = (a", b", c" , d" , e", f") are constants. If we write 
(Y, y, x) as {x^ n+1 \ x^ n \ x^ n ~ 1 '), this is nothing but the symmetric case of the well known 
QRT equation [EEJ CE2] 



(n+l) 



V(x (n) )p"(x (n) ) - p'(x {n) )r]"(x {n) ) -ar( n-1 )(/£/(a;W)0"(a; (n) ) - 0'(^ (n) )p"(^ (n) )) 



0{xW)4J'(x( n )) - , (x( n ))p"(x(™)) - x(«- 1 )^'(iW)i]"(iW) - r/'(x("))0"(a;( n )] 

(49) 

The map fj48|) has an invariant 

(p\x)y 2 + + p\x) 

{,V) (t)"{x)y 2 + i'{x)y + p"{xy 1 ' 

hence it can be reduced to one dimensional map x — > X. The calculation of the 2nd 
elimination ideal is rather trivial in this case. If y and Y are eliminated by using the 
invariant H (x, y) = h, the result we obtain is 

(f)(x)X 2 + rj(x)X + p{x) = (51) 

where 

4>(x) := ax 2 + bx + c, i](x) := bx 2 + (d — 2c)x + e, p{x) := cx 2 + ex + f, 

with 

q = q' + hq. 

If we identify q = (a,b,c,d,e, f) with the one of (1551) . the map (IBTj) is exactly the bi- 
quadratic map (1571) . 



In the theory of the QRT map the formula (15T!) is called an invariant curve. The 
problem of solving the equation fHHl) is now converted to finding the coefficients of (l5Tj) 
iteratively. Our general formula (1471) gives us the explicit expressions of the invariant 
varieties of the symmetric QRT map ( ]48lh 

Recently Tsuda [20] discussed the QRT map from a geometrical viewpoint and ob- 
tained many interesting results. In particular all periodic maps of QRT were classified 
corresponding to particular combinations of the parameters. 

4.4.2 3d Lotka-Volterra map 

After the elimination of x 2 , £3 and X 2 ,X 3 from (}T4"j) by using ffT5]) and denoting x\ = x 
we see that the 4th elimination ideal is given by the biquadratic function ( |38l) with the 
coefficients 

a — r + 1, b = s — 2r — 1, c = r — s, 
d = s 2 + rs + 5r — 2s + 1, e = —r(s + l), f = 0. 

From the general argument we obtain the result of first iteration simply substituting 
these data into (14T1) . They are given by 



= 


(s + r 


) 2 s(r 2 — rs 2 — s — 3rs), 


b {2) = 


( s + r 


) 2 s(2r 2 s + s + 5rs - 2r 2 - r 3 - s 2 ), 


c (2) = 


(s-r] 


)s(r + l)(s 2 — r 2 s — 3rs — r), 


d< 2 > = 


(s + r 


) 2 s(2rs 2 + 2s 2 - 3r 2 s - 8rs - s + r 3 s + r 4 + 5r 3 + 6r 2 - s 3 


e m = 


(s + r 


) 2 rs(s 2 — rs + s — 2r — r 3 — 2r 2 ), 


f(2) = 


(s + r 


) r s(r — rs + r + s + s + 1). 



Notice that all parameters apart from are proportional to a common factor (s + l) 2 . 
This is a result which we cannot derive generally from (f4"Tl) since they are not factorized. 
On the other hand this factor is exactly the one we expect from our previous result ffl6l) . 
although this is a special behaviour of the 3d Lotka-Volterra map. 

The periodicity conditions of period 3 and higher can be derived similarly and the 
associated functions 7®, 7^ and 7^(q) can be read off from ( 1471) directly. We find 

7®(r, s) = r 2 + s 2 — rs + r + s + 1, 

7 ( - 4 - ) (r, s) = 3rs + s + s 3 — 3s 2 r + r 3 s + 6r 2 s — r 3 , 

7 (5) ^ s ) = r 3 s 4 _ r 3 s 2 _ gr 4 s 5 + 1Qr 3 s 6 + 3g 5 r + ^6 + ^5 + 3^4 _ 3r 5 s 3 

_g r 4 s 3 _ r 6 s 3 + 3r 5 s 4 + s 4 + + + ^ g 7 + g 7 + ^j^^ 

-3s e r - r 3 s 5 + 21rV - 10r 3 s 3 - 6rs 7 + s 8 . 

These include again precisely the conditions expected from our direct calculations (flTI) . 
dH]) and (USD. 

It will be worthwhile to mention about the uniqueness of the reduction. If we had 
reduced the map ( |T4l) by eliminating x 2 and X3 one by one in different order we should 
have obtained two different maps. They correspond to the exchange of X^- 1 ' and X^ 1 '. 
This difference becomes irrelevant in our procedure since the 4th elimination ideal (|38|) 
combines them together. 

The 4d Lotka-Volterra map ( |20l) also satisfies d — p= 1. This map, however, does not 
reduce to the biquadratic map, but to a biquartic one. 



4.4.3 Painleve IV map 



The g-Painleve map (I2T!) can be converted into the biquadratic map (|37j) when the pa- 
rameters are fixed at a\ — 012 = c«3 = 1- m fact, after the elimination of X2 and x% we 
find 

a = s + 1, b = — r — 1, c = 0, 
d = r 2 + 4r + sr — s + 1, e = — r 2 — r, / = r 2 — rs, 

which reproduce (|25|) . 

The g-Painleve V map (1261) at «i = 0:2 = a 3 = ol\ = 1 does not reduce to the bi- 
quadratic map but to a biquartic map. 



5 Transition between Integrable and Nonintegrable 
Maps 

We have not discussed, so far, the behaviour of periodic points of nonintegrable maps. The 
purpose of this section is to clarify the transition between integrable and nonintegrable 
maps. In particular we are interested in how a set of isolated periodic points in one side 
turns to an invariant variety in other side. 

It will be useful if there is a simple model which interpolates the integrable and nonin- 
tegrable maps, such that we can study explicitly how the transition of the periodic points 
takes place between the two maps. The method developed in §4.1 provides us many such 
maps, in the sense that the study of one low dimensional map is sufficient to know the 
behavior of many higher dimensional maps asssociated to it, irrespective whether the map 
is integrable or not. 

To be specific we study, in this section, the map (I3T1) which is integrable only in the 
limit c —>■ 0. All information of higher dimensions are included in the parameters a, b, c 
and h. But we assume c is a constant parameter so that we can control the value of c by 
hand without changing the dependence of the higher dimensional coordinates. Moreover 
it is convenient to introduce the new variable z and convert the map ( 13"TT) to the normal 
form of the one dimensional rational maps of degree 2 

2 z = z TTTz' < 52 > 

which has been studied intensively in the literature E] . The correspondence between 
x and the new variable z is given by 

— h + 2hab — hac)q_ — q^_)z + (1 — h + 2hab — hac)q_ + 



x = 

((2/ic — abch — b — bh)q_ + bq\ )z + (2/ic — abch — b — bh)q^ — bq\ ' 

and the new parameters are related by 

(2ch — b)q 2 + — (2hc — abch — b — bh)q_ 

bq\ - (2hc - abch - b - bh)q_ ^ 
, (2ch - b)q\ + (2hc - abch - b - bh)q_ 

bal + (2hc - abch -b- bh)a- ^ ' 



where 

g| = (1 - h - hacf ± Aah(b - he). 

The map ( |52|) is integrable when A A' = 1, hence c = 0, but not integrable otherwise. 
Repeating the map ( 1521) n times the image of z can be written as 

y{n) _ j Z (n-l) ~*~ /-kkN 

A n \A _1 + \A _1 + Z0) J [X- 1 + Z( 2 ) J " ' ^A^ 1 + J ' 1 J 

If is the degree of the rational polynomial Z^ n_1 \ we see from (1551) that k n = 2k n -i, 
hence k n = 2 n . Periodic points are obtained by solving = z. When n is a prime 
number, the number of the periodic points of period n is k n — 2 = 2 n — 2. Here 
subtraction of 2 corresponds to two fixed points. If n is not prime the contribution from 
the divisors must be also subtracted. In this way we find that the number increases as 
fast as 

(#2, #3, #4, • • ■) = (2, 6, 12, 30, 48, 126, 240, 504, ■ ■ ■)• 

In the integrable limit A A' = 1, this increase ceases owing to the cancellation of factors 
in the numerator and the denominator as seen in the second expression (|56j) , and the map 
becomes 

Z (n) = X -n z ( 57 ) 

It is apparent that there is no periodic point unless A satisfies A™ = 1 or 

A = e i 27r fc/n A; = 1,2, ...,n. 



We now recall that A is a function of the invariants of a higher dimensional map. Let us 
write the right hand side of (|53|) as X(h,a,b,c). Then the invariant varieties of periodic 
points of the map are given by 

v [n) = {x | A(/i,a,6,0) n - 1 = 0}, n = 2,3,4,.... (58) 
which provides another expression of (1361) . 



All isolated periodic points suddenly disappear in the integrable limit, and the in- 
variant varieties fl58|) appear. In order to uncover this trick of the transition we must 
know where the isolated periodic points were right before they disappeared. After some 
manipulation we find that the right hand side of (|56l) admits the following expression 

n-2 

I] (^ + A fe ) 2 - fc - 2 + (AA'-l)P n 

y(n) _ k=-l 

\n n-2 

II (* + * k f M + (AA' - l)Q n 
fc=-l 

where P n and Q n are polynomials of z. The periodicity condition Z^ = z requires 



n-2 

(A™ - 1) J] (* + A fc ) 2n " fc - 2 = (AA' - l)(P n - \ n Q n ) 

k=-l 



(59) 



to hold. 

When AA' = 1 we obtain solutions of (1591) at 

z = -l/X,-l,-X,---,-X n - 2 (60) 

if A n 7^ 1. These points are, however, certainly not the periodic points on the invariant 
varieties ( 1581) . Suppose AA' — 1 is very small but not zero. Then the periodicity condition 
f l59|) is satisfied iff x is very close to one of the points satisfying fl60l . Therefore the 
periodic points were in the neighbourhood of these points right before A A' reached 1. 



Finally we consider the connection of our result with the study of the Julia set of the 
map (|52|) . The Julia set can be obtained by calculating iteratively the inverse map of (|52|) 
starting from one of repulsive periodic points. Since the inverse map of (152]) 



\{\Z - A' + ^{XZ + A') 2 + 4(1 - XX') Z) 

(61) 



\{XZ - A' - sj{\Z + A') 2 + 4(1 - XX') Z) 



has two branches, the number of the points increases rapidly and they form a fractal 
structure including many sets of periodic points. 

As the map approaches to the integrable limit A A' — > 1, the inverse map simplifies 
significantly turning to 

r xz 

(62) 

1 -1/A 

and one of the branches is frozen at the point z = —1/A. The degeneracy of many maps to 
this one particular point causes a collapse of the Julia set. We have shown in the papers 
[23j [21] that the Julia set of the map (1521) approaches uniformly to the set lim^oo J n , in 
the limit A A' — > 1, where 

J n = {0,-l/A,-l,-A,---,-A"- 2 }. 

Therefore the points of (1601) are "fossils" of the Julia set after the collapse. We learn, from 
this fact, that the invariant varieties of periodic points are created independent from the 
Julia set in the integrable limit. The Julia set itself ceases to be a set of periodic points 
owing to the degeneracy of the critical points at z — —1/A. 

A simpler model which presents similar property that we discussed in this section 
was invented by S. Onozawa. Using that model he explored more detail of the transition 
between integrable and nonintegrable maps from a different point of view (25]. 



6 Concluding Remarks 

Throughout this paper we have examined many examples of invariant varieties of periodic 
points in order to support our conjecture in §1. To conclude this paper we would like to 
summarize some remarks on our study. 



1. In contrast to the continuous time Hamiltonian flow, the number of invariants does 
not play the critical role to decide integrability of a map. In fact there are higher 
dimensional maps which have many invariants but reduce to nonintegrable maps, 
while some maps, like the q-Painleve maps, have no invariants but are integrable. 

2. If a periodic point of some period is isolated, the periodic map can start only from 
certain limited points. When the point is on an invariant variety of periodic points, 
on the other hand, the map can start from arbitrary point of the variety. Therefore 
the conjecture in §1 can be rephrased as "an existence of movable periodic points, 
dependent on initial values, indicates integrability of the map." 

3. The Poincare sections of a periodic orbit of an integrable system produce a smooth 
curve in a plane of the phase space as initial values are changed continuously. If we 
consider the intersections as a sequence of a discrete map, we can think of the curve 
as an invariant variety of periodic points. When the map is perturbed, our lemma 
requires that all points on the invariant variety are frozen by the Birkhoff's fixed 
point theorem. 

4. It is well known that some higher dimensional nonintegrable maps are characterized 
by strange attractors in which a dense set of periodic points of all periods are 
condenced. As we consider the attractor in the complex space, it is embedded in 
the Julia set in some cases [9], [26] . Therefore it is a different object from the invariant 
variety of periodic points, since the latter is a variety formed by a set of periodic 
points of a single period. 

5. We have concerned in this paper with a sufficient condition for the integrability 
of a rational map. The necessary condition for the integrability will be given iff 
the precise notion of integrability is known. Our theorem can not discriminate an 
integrable map from nonintegrable one if the periodicity conditions are uncorrelated 
in all periods. The study of the g-Painleve maps in §3-3 shows the difficulty of 
clarifying the border of these two regimes. 

6. In order to clarify the phenomenon which takes place at the boder between inte- 
garable and nonintegrable regimes we have investigated in detail the transition of 
a map from one side to another in §5. It is found that the source of the periodic 
points are clearly different. Namely the periodicity conditions are imposed on the 
variables to determine their position if the map is nonintegable, while, if the map is 
integrable, they are imposed on the invariants to determine their relations among 
themselves. This argument, however, does not apply to an integrable map which has 
no invariant. We also would like to mention that there are some subtleties involved 
in defining integrability even when all periodic points are neutral [271 128] . 

Acknowledgements 

The idea of this study arose from discussions with Dr. Katsuhiko Yoshida some years 
ago. We would like to thank him and also Mr. Show Onozawa for many interesting discus- 
sions. We would like to express our special thanks to Prof. Martin Guest for suggesions 
concerning the exposition. We also thank Prof. Akira Shudo for many valuable comments 
especially on the topics related to his own works. 



Appendix: Invariants of the Lotka-Volterra Maps 



Let us define the matrices 



R(t) 



L(t) 



( 1 - X? 






1 

/ 

X® 




1 - X 



(') 






1 





1 

1 -X, 


X. 



(t) 

d-l 





1 



l-Xf J 



(0 



x 



it) 



V 









1 

X. 







1 \ 

(*) 



it) 

d-l 



A straightforward calculation will show that the Lotka-Volterra equations (|T2l are equiv- 
alent to the following matrix formula, 

L{t + l)R{t + 1) = R(t)L(t). 

To find invariants we define 

••• 

1 
p 2 1 1 ••• 
p 3 1 1 ••• 



/ 1 

1 



A(t) := L{t)R{t) 



Pd 1 \ 
pi 






V 



1 1 

pd_i 1 1 / 



where we used the notation j?j := X) (1 — Since 

A(t + 1) = L(t + l)R(t + 1) = i?(t)A(t)i? _1 (t), 
eigenvalues of A are invariant. If we write 

det(A-A) = (-l) d - 1 X:^(A-l) fc ! 

fc=0 

the set of coefficients H , Hi, H d is also invariant. Comparing both sides we find 

•l) d PiP2 • • -Pd, k = 

■ Pjk , k = l,2,...,[d/2] 

0, k=[d/2] + l,...,d-l 

— 1, k = d. 



Hh 



1 



Here [d/2] = d/2 if d is even and [cZ/2] = (d— 1)/2 if d is odd. The prime in the summation 
2' means that the summation must be taken over all possible combinations ji,j2i ■■■ijk 
but excluding direct neighbours. Since H can be represented by other H^s and 



(*) y(t) 

2 



it is convenient to use r instead of H 



■X. 



it) 

d i 
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